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Abstract. We consider the applications of functional renormalisation group to few and many-body systems. As 
an application to the few-body dynamics we study the ratio between the fermion-fermion scattering length and 
the dimer-dimer scattering length for systems of few nonrelativistic fermions. We find a strong dependence on the 
cutoff function used in the renormalisation flow for a two-body truncation of the action. Adding a simple three- 
body term substantially reduces this dependence. In the context of many-body physics we study the dynamics 
of both symmetric and asymmetric many-fermion systems using the same functional renormalisation technique. 
It is demonstrated that functional renormalisation group gives sensible and reliable results and provides a solid 
theoretical ground for the future studies. Open questions as well as lines of further developments are discussed. 



Nonperturbative/Exact Renormalisation Group (ERG) 
has become a popular tool to study strongly interacting dy- 
namics IH. It has been used to study a variety of physical 
systems, from nonrelativistic two-body |2| and many -body 
O systems to effective quark models [4J and gauge theo- 
ries 0. The main object of the ERG is the scale-dependent 
average effective action r k , where k is an auxiliary running 
scale (see the reviews O for detailed discussion). The evo- 
lution of A with the running scale is described by a flow 
equation. The average effective action at scale k 2 contains 
the fluctuations with momenta q 2 larger than k 2 . Fluctua- 
tions with momenta smaller than k 2 are integrated out. In 
the limit k — > all the fluctuations are included and the 
full effective action is recovered. To integrate out modes 
with q 2 < k 2 one introduces a cutoff function Rk(q) in path 
integral. This function should vanish in the limit k — > to 
ensure the physical results to be cutoff independent and be- 
have like k a , a > when q — > 0. In the limit of large k the 
average effective action should turn into the basic action. 

In this paper we consider several applications of ERG 
to the physical systems. Namely, we calculate the scatter- 
ing length of two dimers - bound states of two fermions, 
which play an important role in dynamics of ultra-cold and 
dilute Fermi gases. We also consider both symmetric and 
asymmetric many-fermion systems with pairing interac- 
tion and analyse some physical characteristic of these sys- 
tems such as pairing gap, chemical potential and ground 
state energy. 

The field of ultra-cold Fermi gases has become one of 
the hot topics in atomic physics. An important feature of 
such systems is superfluidity which is the result of attrac- 
tive fermion-fermion interactions leading to pairing. Re- 
cent advances in using Feshbach resonances allow a tuning 
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of the fermion-fermion scattering length a. For negative 
scattering length we get the small-coupling BCS state. For 
positive values of a bound states of two fermions (dimers) 
form and these can be arranged into a Bose - Einstein con- 
densate (BEC) |7 1 ■ The size of dimers is determined by the 
fermion-fermion scattering length and their binding energy 
is of order 1/a 2 . 

The exact relation between dimer-dimer and fermion- 
fermion scattering lengths a# — 0.6a was established in 
1 8 ] by solving the Schrodinger equation for two composite 
bosons interacting with an attractive zero-range potential. 
This method, however, is difficult to extend to the many- 
body case. Therefore, it is useful to study the ratio ag/a 
in an approach which can be used both for few and many- 
body problems. 

For an exact solution of the functional RG equation we 
expect cut-off independent results for ag- In practice, how- 
ever some dependence of the results on the cut-off is in- 
evitable. We can use this dependence as a measure of the 
completeness of our parametrisation. We will see that the 
ansatz used for the effective action with two-body interac- 
tion is too crude and needs to be modified by including the 
term with the three-body interaction. 

The flow of the scale-dependent effective action satis- 
fies 

d k r = -LsTr[(d k R)(r™~Rr 1 ]. (1) 

where 7" (2) is the second functional derivative taken with 
respect to the fields entering the action, and R is an oper- 
ator that drives the RG evolution. The operation STr de- 
notes the supertrace [9| taken over both energy-momentum 
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variables and internal indices and is defined as 

S Tr A a b J f J = Tr(A BB ) - Tr(A FF ). (2) 

In the case of mixed boson-fermion systems d^Riq) acts 
on both boson or fermion degrees of freedom. As men- 
tioned above the flow equation is a functional differential 
equation, and since there are no general methods to solve 
such equations numerically, we must resort to approxima- 
tions. One approach is to parametrise the effective action 
by a finite set of coefficients, turning the evolution into 
a system of coupled ordinary differential equations which 
can then be solved numerically. In this paper, where we 
study possible truncations for the system of a few inter- 
acting s-wave fermions, our choice of ansatz for effective 
action is motivated by both studies of many-body systems 
and few-fermion systems in the effective field theory 
(EFT) approach [ 10 1 . The study of such few-body physics 
in the framework of ERG has just begun. The famous Sko- 
rnyakov and Ter-Martirosyan equation has been rederived 
in Ref. ifTTl and Efimov type of physics has been consid- 
ered in Ref. IE) (see also EH). Here we follow Ref.lfl4l 

We first consider the case when fermions interact pair- 
wise. This case was previously considered by the Heidel- 
berg group [ 15 1 and more recently in 0, where many the 
technical details can be found. Here we just give a brief ac- 
count of the formalism, concentrating on the dimer-dimer 
scattering length. We reproduce the results from Ref. lfl5l 
which allow us to establish the framework for examin- 
ing the cutoff scale dependence of the results. At some 
starting scale we demand the microscopic action to be a 
purely fermionic theory with the contact 4-fermion interac- 
tion without derivatives. This kind of interaction has been 
extensively used in the EFT-based studies of nucleon-nucleon 
forces iflOl . It is convenient to modify the theory by in- 
troducing an auxiliary composite boson field. Using the 
Hubbard-Stratonovich transformation the 4-fermion inter- 
action gets replaced by a "Yukawa" type of coupling be- 
tween the fermions and an auxiliary boson. A kinetic term 
for the boson is then generated in the RG evolution. We 
also add a local three-body interaction, which we expect 
to be generated during the evolution, to the standard two- 
body result. This can be captured in a following parametri- 
sation of the effective action, 



r[<fr, r , 4>, <t>\ k] = r 2b [^, if/ T , 0, <f>\ k] 

-A j t/ 4 X(/r t (x)^ t (x)0(x)(/r(x), (3) 



where 



j (fx j d 4 x' 4>\x)n(x, x'; k)4>(x) 
1 



+r(x)\id t + — V z )Hx) 



-g\ - ^{x)cr 2 ifr{xW{x) - - f (x)cr 2 f [ (x)<p(x) 



4) 



We shall first concentrate on the two-body part F 2 b- 
Here IJ(x, x';&) is the scale-dependent boson self-energy. 
The evolution of this self-energy is given by 



d k n(x,x';k) = 



S(f>(x')Sfl(x) 



dkr\<p=o, 



(5) 



but from here-on we shall express all evolution in momen- 
tum space. Note that only fermion loops, which only de- 
pend on the fermion cut-off Rf and its derivatives, con- 
tribute to the evolution of the boson self-energy in vacuum. 
Integrating the ERG equation with the fermionic cut-off 



k 2 -a 2 



(6) 



and the renormalisation condition that the constant term 
in I7(P(), P; 0) reproduces the inverse of the zero energy T 
matrix, we find 

g 2 M [4 n 16 / P 2 \ P 3 

n(P ,P;K) = — n K+- + — \MP r 

v ; 4n 2 [ 3 a 3K\ 2 J 24K 2 

(7) 

The on-shell fermion-fermion scattering amplitude in the 
physical limit k — > is given by 



-L = ln(P ,P;oi 
Tip) g 2 



(8) 



where p = y/2MP - P 2 /2 is the relative momentum of 
two fermions and Pq{P) denote the total energy (momen- 
tum) flowing through the system. From the gradient expan- 
sion, we define boson wave-function and mass renormali- 
sation factors by 



Z,(*) = — 77(P ,P;*) 
oPo 



(9) 



p =e D ,p=o 



and 



AM 



Z m (k) 



dP 2 



n(P ,P;k) 



(10) 



p =e d ,p=o 



Here &d = -I /(Ma 2 ) which is the bound-state energy of 
a pair of fermions. Note that these renormalisation factors 
are only identical in vacuum for a limited subset of cutoff 
functions preserving Galilean invariance to otherwise the 
identity Z^(k) = Z m (k) holds only in the physical limit k — > 
and the evolution of these renormalisation factors should 
be calculated separately. 

The evolution of the boson-boson scattering amplitude 
follows from 



(In) 



■d k u 2 (6 D , k) = 



S 4 



<5^(fi D ,O)(50t2( 6Dj o) 



(11) 

The evolution can be separated into fermionic and bosonic 
contributions containing OuRf and 8i<Rb, respectively. We 
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first look at the mean-field result, where boson contribu- 
tions are neglected. The evolution of u 2 is then given by 

kUl 4 J (2^[(£F fi (q,Q-£ D /2)] 4 ' 
where E FR (q,k) - th q 1 + Rp{q,k)- Explicit calculations 



give 



"2(0) = — MV a 3 
Ion 



(13) 



where we have again used the sharp cut-off function of the 
form Eq. ©. 

The mean field scattering amplitude at threshold is 

8tt 2u 2 (0) ina 

Trr = a R = r — = . (14) 

BB 2M B Z 2 M ' 



This is the well-known mean-field result a B — 2a |fl6ll 
which is far from the exact value a B — 0.6a [8 J. This im- 
plies that beyond-mean-field effects such as dimer-dimer 
rescattering are important to be considered. To include these 
effects one needs to take into account the boson loops. Af- 
ter some algebra we get 



d k u 2 = d k u 2 \ F + d k u 2 \ B 



Here 



d k u 2 \, 



u\{k) 



s 



d 3 q 



d k R L 



2Zj(k) J (2tt) 3 [E BR (q,k)-6 D y 



and 



„ . 3/ r d\ d k R f 

d ^ = - — J j2nP[E FR ( q ,k)- 



k)-S D /2]' 



(15) 



(16) 



(17) 



where 



and 



E B R(q, k) = — q 1 + — — + - — (18) 



4M Z#(k) Z#(k) 



ui(k) = -n(S D ,0;k). 



(19) 



The boson cut-off function is to be chosen as close as pos- 
sible to the fermionic one, 

R B (q, k) = Z (CBk 2 M q2 () ^ k ~ ?). (2°) 

apart from the addition of a parameter c B , which sets the 
relative scale of the fermion and boson regulators, and a 
factor of Z^. The latter has two main advantages: Firstly, 
it leads to universality, where all contributions to a sin- 
gle evolution equation decay with the same power of k for 
large k and secondly we get a F scaling, where all terms in 
a single evolution equation have the same dependence on 
a F . 

The mean-field result is recovered when c B = 0. The 
limit of c B — > oo leads to a B — > 0. Using c B = 1 gives 
the value of the ratio a B /a = 1.13. Choosing c B — V2 
results in a B /a F - 0.75 as in Ref. lfT5l . In general, the 



results show rather strong dependence on the boson scale 
parameter c B (see Fig. [TJ. This is unwanted as it means 
that the physical results depend strongly on the choice of 
regulator. As announced above we shall now consider the 
effect of the three-body force (??). 

The equations for u \ and remain unchanged but the 
one for u 2 gets modified by the term describing the 3-body 
interaction : 



d k u 2 =-2Ag J — 



d k R 



k*F 



f [E FR (q,k)-S D /2f 



(21) 



where we denote E FR ,p B = E FR -P B /2 and E BR p B = E BR — 
P B . The evolution equation for A is defined by an expansion 
at the energy of the bound state pole for bosons, and half 
of that for fermions, 



cM = --- 



5 A STr[d k R(r {2) - R)- x \ 



2 5<pH6 D , 0)S<f>(S D , 0)6ifr1(6 D /2, 0)5ifr(P B /2, 0) 

(22) 

There are three distinct contributions to the running of 
A coming from a ladder, triangle and box diagrams. We 
denote the corresponding driving terms as D/,, D\ and D t , 
splitting the last two terms in fermionic and bosonic con- 
tributions. After evaluation of traces and contour integrals 
we get 



2 P d\ 
J (2tt 



d k R F Z<t, + d k R B 



q 

) 3 (E FRyPls Z0 + E B 



D 



Df - g 2 A 



r d\ 

J (2^ 

r d\ 

J Q*. 



d 3 q d k R F {E BR p B + 2Z <p E FRiPb ) 



) E FR p B (.E FR p B Z/p + E br?Pb ) 2 



d k R B 



nF _9* f d s q 
h ~4j (2n? 

B _g 4 f d\ 



) 3 Efr,p b (Efr,p b Z<p + E BR ,p B ) 2 



d 3 q d k R F (2E BRt p B + 2>Z$E FR f B ) 



E FR j, b (Efr,p b Zq + E br , Pb ) 2 



d k R B 



Efr,p b (Efr,p b Z,i> + E br , Pb ) 2 



(23) 



(24) 



(25) 



(26) 



(27) 



We denote E fr , Pb = E FR - P B /2 and E B r,p b = E BR - P B 
with P B being some external energy. The evolution equa- 
tion for A is defined by an expansion at the energy of the 
bound state pole for bosons, and half of that for fermions, 
Note that the the evolution of the three-fermion term A does 
not depend on the four-fermion term interaction u 2 , but u 2 
depends on A. We now need an initial condition for A. At 
infinite k, there is no fundamental three-fermion interac- 
tion, and thus A = 0. To simplify the calculations we have 
assumed that A = is also zero at the starting scale, since 
A behaves like l/k 2 this can be used reliably. 

Now we turn to the results. We first note that the results 
are numerically independent of the choice of starting scale 
provided it is chosen to be at least ka F ^ 100. For k » 
1 /a the system is in the "scaling regime" and evolves near 
the fixed point until the scale becomes comparable with 
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Fig. 1. Ratio of boson-boson to fermion-fermion scattering length 
as a function of the scale parameter c B . The blue (solid) curve 
shows results without 3 body forces, the red (dashed) curve con- 
tains the effect of a 3-body force. 



1 /a. This can easily be seen if we switch to dimensionless 
coupling k 2 A and analyse its evolution. 

We have found that the ratio of as I a decreases when 
the three-body term is included. For example, assuming 
eg — 1 leads to the value ag/a = 0.74 compared with 
ag I a = 1.13 without inclusion of the the three-body term, 
and the choice eg = V2 as used in [ 15 1 gives a^ja = 0.69. 
We have shown in Fig. 1 the behaviour of ag/a as a function 
of the scale c#. 

The upper curve corresponds to the calculations with- 
out three-body term and the lower curve gives the results 
of the full calculations with three-body term included. It 
is seen from Fig.l that taking into account the three-body 
forces results in significantly weaker dependence of as I a 
on the scale eg. Note that at large eg, which corresponds 
to integrating out fermions first, the two curves approach 
each other. Indeed, in this limit the dominant contribution 
to cis/a comes from the bosonic part (proportional to Rb) 
of the equation for m 2 which does not depend on the three- 
body coupling A. On the other hand, with small cb the main 
contribution comes from the fermion loops which get mod- 
ified when the three-body coupling A is included. We ex- 
pect this picture to be qualitatively correct for any type of 
boson regulator although the quantitative details will vary 
depending on the functional form of the regulator used. 

In spite of showing a clear improvement over the cal- 
culations with the two-body interactions only, adding the 
simplest possible three-body term is not enough to ensure 
the complete (or almost complete) independence of the re- 
sults from the scale parameter eg in the region eg — 1. It 
is worth emphasising again that, as long as any truncation 
of the effective action is made, the results will never be 
completely independent of the scale c#. Instead, we expect 
some "stability window" where the slope of as/a is small. 
It is somewhat similar to the QCD sum rules JT7), where 
the physical observables are expected to be independent of 
an auxiliary parameter called Borel mass. It seems likely 
that further extensions of the effective action will result in 
stability of the results with respect to the variations of c# in 
wider region. Such extensions may include both four-body 
interactions as well as energy and/or momentum depen- 



dent three-body forces. Work along these lines is now in 
progress. 

We next consider the applications of the ERG to the 
finite density many-fermion systems. It is well established 
that in the medium even a weak attraction between fermions 
may lead to the intrinsically nonperturbative phenomenon, 
the superfluidity, characterised by rearrangement of the ground 
state and appearance of the gap in the spectrum. The fermions 
form correlated pairs which, depending on the strength of 
the interaction, may lead to different physical regimes. The 
weak coupling regime (BCS phase) corresponds to a pair 
with the spatial size much larger then the radius of the in- 
teraction so that no actual bound two-body subsystem is 
formed, while in the strong regime corresponding to the 
Bose-Einstein Condensation (BEC) the fermion pairs form 
compact deeply bound two-body states. We demand that at 
high scale our theory be a purely fermionic theory with the 
contact interaction described by the lagrangian 

■Ci = -\ Co (iAV 2 iA tT ) (<aV 2 (A) , (28) 

which is identical to one, used to calculate the fermion- 
fermion scattering in vacuum so we can link the vacuum 
and in-medium ERG calculations. This is what makes the 
approach universal. Since we are interested in the appear- 
ance of the correlated fermion pairs in a physical ground 
state, we need to parametrise our effective action in terms 
of corresponding variables. A natural way to do this is to 
introduce a boson field whose vacuum expectation value 
(VEV) describes this correlated pair and study the evolu- 
tion of this effective degrees of freedom. At the start of the 
RG evolution, the boson field is not dynamical and is in- 
troduced through a Hubbard-Stratonovich transformation 
of the four-point interaction. As we integrate out more and 
more of the fermion degrees of freedom by running k to 
lower values the dynamical term in the bosonic effective 
action is generated. 

The corresponding ansatz for the boson-fermion effec- 
tive action can be written as 

v(z$(id t +2fi) + ^ 2 y-U(<p 



-g <A T o- 2 iAf - l - ^V 2 ^ tT 0J 



Here M is the mass of the fermions in vacuum and the 
factor l/2m in the boson kinetic term is chosen simply to 
make Z m dimensionless. The couplings, the chemical po- 
tential the wave-function renormalisations and the 
kinetic-mass renormalisations Z,„,m all run with k, the scale 
of the regulator. The bosons are , in principle coupled to the 
chemical potential via a quadratic term in (f>, but this can be 
absorbed into the potential by defining U = U - IfjZ^'tp. 
We expand this potential about its minimum, <p^<p = po, so 
that the coefficients m, are defined at p — po, 



1 

U(p) = u + u l (p-p )+- u 2 (p-p Q Y 



1 , 

+ -M 3 (p-p ) +■ 

o 



(30) 
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where we have introduced p = <ffi<f>. The phase of the sys- 
tem is determined by the coefficient U\. In the symmetric 
phase we have po — so that the expansion takes the form 



0.15 



1 2 

U(p) = u + uip+ - u 2 p + 



(31) 



The potential in the condensed phase can be simplified to 



1 , 
U(p) = u + - u 2 (p - po) + 



(32) 



In our current work we shall truncate this potential at quar- 
tic order in the field (order p 2 ). We treat the wave function 
renormalisation factor for the bosons in the same way, ex- 
panding it about p = po as 



Z^(p) = z^o + z*i (fi - Po) + 



(33) 



The other couplings and renormalisation factors can be 
treated similarly. 

The fermions are not dressed at this point, the bosons 
are just auxiliary fields and therefore we can assume that 
Z#(K) = 1 and Zm{K) = 1 . After rather lengthy algebra 
one can obtain the explicit expressions for the couplings 
in the effective action. As an example we show below the 
corresponding ERG equation for the effective potential. 



d k U = 



1 C d 3 q 



4 



; dkRf 



2Z 



i r d\ 



E FR +A* 



-'BR 



4 



: d k R B . (34) 



^ BR ^B 



where 



E B R(q,k) = - ! ^-q 2 +u 1 +u 2 (2f 
2m 



-po)+R B (q, k), V B = u 2 <f> 
(35) 

The evolution equations for the coefficients in our expan- 
sion of the effective potential are obtained from the deriva- 
tives of dkU with respect to p = <p^(f>. 

The initial conditions for them can be fixed by demand- 
ing that at high scale the evolution of the system is close to 
the vacuum one. One could introduce a cut-off function that 
tends to a p^-independent form for K » p F . However in 
practice a modification of the renormalisation procedure is 
more convenient. In the region K » p F , we can ignore bo- 
son loops. The evolution of quantities such as u\(p F ,K), 
u 2 (pf,K), Z${p F ,K) and Z„,{p F ,K) is thus similar to the 
vacuum case, except for the different cut-off. This allows 
us to define u\(p F , K) to be 



ui(p F ,K) 
9 2 



m i r d 3 i 

~ ~4^ + 2 J (2^r 




k (fm ) 



Fig. 2. Numerical solution of the evolution equations for infinite 
scattering length.The thinner curve (orange online) corresponds 
to the mean-field calculations and the thicker one (blue online) 
includes boson loops. The difference however is noticeable only 
for the evolution of u 2 (graph in the middle). 



Rf (q, 0, k) for k pp . It ensures that our renormalised pa- 
rameter u\(pf,K), defined using Rf(q,pf,k) for large k, 
corresponds to the physical scattering length in vacuum. 

The initial values for u 2 (pf, K), Z^,(pp, K) and Z„,(pf, K) 
can be determined using similar procedures, although this 
is not so crucial since these quantities do not contain lin- 
early divergent pieces and so all their pp -dependence is 
suppressed by powers of Pf/K. One convenient choice is 
to take their starting values to be zero at some large but fi- 
nite scale K. An alternative is to require that they tend to 
zero as K — > oo. 

The initial condition for the energy density is most con- 
veniently expressed in terms of So which, in the symmetric 
phase, is simply given by the energy of a free Fermi gas, 
measured relative to the chemical potential, and so its ini- 
tial value is just 



u (K) = 2 



r a\ 

J (2n) : 



E FR (q,p F ,0)6(p F -q). 



(37) 



1 



sgn(q - p F ) 



E FR (q, 0,0) E FR (q,p F ,K) 
(36) 

This expression can be thought of as being generated by 
the vacuum evolution using a modified cut-off that inter- 
polates smoothly between R F (q, p F ,k) for k » p F and 



Now we turn to the results which are shown on Fig. 
2. First we note that the value of the physical gap is prac- 
tically independent of the starting scale K provided K > 

5/m" 1 . 

At starting scale the system is in the symmetric phase 
and remains in this phase until u \ hits zero at k crlt 1 .2fm~ l 
where the artificial second order phase transition to a bro- 
ken phase occurs and the energy gap is formed. Already 
at k =i 0.5 the running scale has essentially no effect on 
the gap. The results obtained for the gap correspond to the 
case of the infinite negative scattering length. To study the 
B'CS-BEC crossover we have to solve the evolution equa- 
tion for a wider range of the scattering lengths, including 
the positive values. The corresponding results for the evo- 
lution of chemical potential as a function of the parameter 
p F a are shown in Fig. 3. 
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Fig. 3. Evolution of the chemical potential 

While the vacuum scattering length is large and neg- 
ative, the system is in the BCS-like phase with positive 
chemical potential, whereas if the scattering length is cho- 
sen to be large and positive reflecting the existence of a 
bound state near threshold the system ends up being the 
collection of weakly overlapping tightly bound pairs with 
negative chemical potential. 

As in the case of dimer-dimer scattering we expect here 
some dependence on the regulator used. By analogy with 
the vacuum case the inclusion of the three-body forces will 
probably make the results more stable although the corre- 
sponding calculations are much more involved compared 
to the vacuum case. Work along this line is in progress as 
well. 

We have also calculated the ground state energy of the 
many-fermion system in the unitary regime - the idealised 
case of the infinite fermionic scattering length ifTSl . In this 
regime the dynamics of the many-body system becomes 
independent of the microscopical details of the underlying 
interaction between two fermions. In the a — > -oo limit 
the ground state energy per particle is proportional to that 
of the non-interacting Fermi gas. 

E GS =tE FG =Z^=^E F , (38) 

where M and k F are the fermion mass and Fermi momen- 
tum correspondingly and £ is the universal proportionality 
constant, which does not depend on the details of the in- 
teraction or type of fermions. The other dimensional char- 
acteristics of the cold Fermi-gas in the UR such as paring 
energy A or chemical potential // can also be represented in 
the same way 

yu = 7]E FC , A = eE FG (39) 

The infinite scattering length implies nonperturbative treat- 
ment that why the ERG formalism is a very useful tool 
here. We found the values of 0.62,0.96 and 1.11 for the 
universal coefficients ^, r\ and e correspondingly. 

The obtained value for £ is close to the experimen- 
tal data from [19], £ = 0.74(7). The other measurements 



give f = 0.34(15) and rj = 0.6(15)03; f = 0.32^ 

and rj = 0.53+°'Jo EQ; £ = 0.46(5) and rj = 0.77(5)ES; 
£ = 0.51(4) and rj = 0.85(4) (23]]. The results of theoretical 
calculations are summarised in Table I. 



Table 1. Universal coefficients 



Ref 




e 


1 


|24 1 


0.42 


0.9 


0.71 


m 


0.22 






ED 


0.41 




0.7 


(26) 


0.44 


0.93 




(28) 




1.03 




(29) 


0.39 






ED 


0.3 


0.66 




EO 


0.55 







As one can see both experiment and numerical simu- 
lations do not provide the coherent value of the £ constant 
so it is difficult to judge the quality of the numerical es- 
timates provided by the ERG calculations. One may only 
conclude that the ERG approach leads to the sensible val- 
ues of the universal coefficients consistent with the exper- 
iment and lattice calculation but more detailed compari- 
son can be done when more accurate date are obtained. We 
note, however that the value of the universal parameters are 
still somewhat higher then the "world average". One pos- 
sible cause could be the neglection of the screening effects 
(32) which are known to decrease the values of the gap 
and energy density. Naive extrapolation of our results us- 
ing the known value of the Gorkov - Melik-Barkhudarov's 
correction |32| indeed brings the values of the universal 
coefficients closer to the "world average" of the lattice and 
experimental data. Clearly, this point requires further anal- 
ysis. 

We now turn to the case of asymmetric many-fermion 
system. This asymmetry can be provided by unequal masses, 
different densities and/or chemical potentials. Understand- 
ing the pairing mechanism in such settings would be of 
immense value for different many fermion systems from 
atomic physics to strongly interacting quark matter. The 
important theoretical issue to be resolved here is the nature 
of the ground state. Several competing states have been 
proposed so far. These include: LOFF (3311 phase, breached- 
pair (BP) superfluidity [34] (or Sarma phase) and mixed 
phase (35). Establishing the true ground state is still an 
open question. It was shown, for example, that LOFF and 
mixed phases are more stable then the Sarma phase in the 
systems of fermions with the mismatched Fermi surfaces 
and with both equal and different masses (36). We concen- 
trate on the case of two fermion species. The ansatz for the 
effective action we use to run the ERG evolution is a sim- 
ple generalisation of the one used for the many-fermion 
system with one type of fermions and consists of summing 
over two types of fermions. Calculating the second func- 
tional derivatives, taking the matrix trace and carrying out 
the pole integration in the loop integrals we get the evolu- 



19 International IUPAP Conference on Few-Body Problems in Physics 



tion equation for U at constant chemical potentials 

i r d 3 a Ef? 

dkU = -~ | ; — [dkR F ,a + d k R F , h ] 



1 r d 3 q E FS 

= ~T 7T^3 I [OkRF,a + d k R F j_ 

2 J W Je%~^7 2 
i r d\ 

2Zl] (2n) 



(40) 



It is worth mentioning that poles in the fermion propa- 
gator occur at 



?o 2 = ~ E * ± 



2 +A 2 . 



(41) 



Here 



Es - {E FR , a +E FRyb )/2, E A - (E FR ^ a -E FRyb )/2, (42) 
and 

E FRJ (q, p^i, k) = q 2 - m + R F (q, k), A 2 = g 2 $<f>. 

(43) 

and we have introduced p^i — J2MifTi, the Fermi momen- 
tum corresponding to the (running) value of At k = 
(R F = 0) in the condensed phase, these become exactly the 
dispersion relations obtained in ||34l where the possibility 
of having the gapless excitations has been discussed. The 
ordinary BCS spectrum can easily be recovered when the 
asymmetry of the system is vanishing (Ea — > 0). The first 
term in the evolution equation for the effective potential de- 
scribes the evolution of the system related to the fermionic 
degrees of freedom whereas the second one takes into ac- 
count the bosonic contribution. The mean field results can 
be recovered if the second term is omitted. In this case the 
equation for the effective potential can be integrated ana- 
lytically and then, by differentiating the effective potential 
with respect to p and setting the derivative equal to zero, 
one can derive the standard mean-field gap equation. 

In this paper we consider the simplified case of two 
fermion species with the different masses and the same 
Fermi momenta 071 . It implies that the chemical poten- 
tials are different. In this situation neither the Sarma phase 
nor the LOFF phase exists and the system experiences the 
BCS pairing depending however on the mass asymmetry. 
The general case of the mismatched Fermi surfaces will 
be discussed in the subsequent publication. For simplicity 
we consider the case of the hypothetical "nuclear" matter 
with short range attractive interaction between two types 
of fermions, light and heavy, and study the behaviour of 
the energy gap as the function of the mass asymmetry. 
We choose the Fermi momentum to be p F = \31fm~ 1 . 
One notes that the formalism is applicable to any type of 
a many-body system with two fermion species from quark 
matter to fermionic atoms so that the hypothetical asym- 
metrical "nuclear" matter is simply chosen as a study case. 
We assume that M a < Mb, where M a is always the mass of 
the physical nucleon. 

First we consider the case of the unitary limit where the 
scattering length a = -oo. The results of our calculations 
for the gap are shown on Fig. 4. 
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Fig. 4. Evolution of the gap in the MF approach (dashed curve) 
and with boson loops (solid curve) in the unitary regime a = -oo 
as a function of a mass asymmetry. 



We see from this figure that increasing mass asymme- 
try leads to a decreasing gap that seems to be a natural re- 
sult. However, the effect of the boson loops is found to be 
small. We found essentially no effect in symmetric phase, 
2 - 4% corrections for the value of the gap in the broken 
phase and even smaller corrections for the chemical poten- 
tial so that one can conclude that the MF approach indeed 
provides the reliable description in the unitary limit for 
both small and large mass asymmetries. It is worth men- 
tioning that the boson contributions are more important for 
the evolution of U2 where they drive U2 to zero as k — > 
making the effective potential convex in agreement with 
the general expectations. This tendency retains in the uni- 
tary regime regardless of the mass asymmetry. 

We have also considered the behaviour of the gap as 
the function of the parameter p F a for the cases of the zero 
asymmetry M a = Mb and the maximal asymmetry Mb = 
10M„. The results are shown on Fig.5. 

One can see from Fig.5 that in the case of zero (or 
small) asymmetry the corrections stemming from boson 
loops are small at all values of the parameter p F a consid- 
ered here (down to p F a = 0.94). On the contrary, when 
Mb = 10M a these corrections, being rather small at p F a > 
2 becomes significant (~ 40%) when the value of p F a de- 
creases down to p F a ~ 1. We found that at p F a ~ 1 the 
effect of boson fluctuations becomes ~ 10% already for 
Mb = 5M a . One can therefore conclude that the regime 
of large mass asymmetries, which starts approximately at 
Mb > 5M a , with moderate scattering length and/or the 
Fermi momenta is the one where the MF description be- 
comes less accurate so that the calculations going beyond 
the MFA are needed. One might expect that the devia- 
tion from the mean field results could even be stronger in 
a general case of a large mass asymmetry and the mis- 
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Fig. 5. Evolution of the gap as a function of the parameter pf a. 
The upper pair of the curves corresponds to the calculations with 
no asymmetry in the MF approach (dashed curve) and with boson 
loops (solid curve) and the lower pair of the curves describes the 
results of calculations with the maximal asymmetry when Mt = 
10 M a 



matched Fermi surfaces but the detailed conclusion can 
only be drawn after the actual calculations are performed. 

In summary, we have performed the ERG analysis of 
a number of fermionic systems. In particular we have con- 
sidered boson-boson scattering length, where bosons are 
just bound states of two fermions. Our study has revealed 
that, starting from rather simple ansatz with two-body in- 
teraction term one can get the results, close to the exact 
value. However, the results turn out quite sensitive to the 
value of cb, the boson scale parameter. We have shown that 
the inclusion of the three-body interaction brings the ratio 
as/a closer to the exact value and significantly reduce the 
sensitivity of the results to the boson scale parameter. We 
have also studied the many-fermion systems with attrac- 
tion leading to pairing phenomena and calculated a number 
the many-body parameters such as gap and its dependence 
on mass asymmetry, chemical potential both in BCS and 
BEC regions and the ground state energy. The directions 
for further improvement of the approach are outlined. 
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